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Abstract 

Symmetries of a differential equations is one of the most important 
concepts in theory of differential equations and physics. One of the 
most prominent equations is KdV (Kortwege-de Vries) equation with 
application in shallow water theory. In this paper we are going to 
explain a particular method for finding symmetries of KdV equation, 
which is called Harrison method. Our tools in this method are Lie 
derivatives and differential forms, which will be discussed in the first 
section more precisely. In second chapter we will have some analysis 
on the solutions of KdV equation and we give a method, which is 
called first integral method for finding the solutions of KdV equation. 
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Introduction 

There are diverse methods for finding symmetries of differential equations. 
One of the most important one is Lie method. In this method we need a lot 
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of basis, such as Lie group theory, prolongation and which are foundations 
of Lie's theory of symmetry groups of differential equations. This is a firm 
method for finding symmetries and it has a lot of applications in differential 
equations theory. Harrison's method is an another one, which does not need 
the above necessaries. Here we act this method on the KdV equation, a very 
famous differential equation, and testing the results, specially with the Lie 
method see [4], [14], [15], [16]. 

1 Harrison's Method 

The method proceeds as follows. We consider a set of partial differential 
equations, defined on a differentiable manifold E with p independent vari- 
ables and q dependent variables, here we assume p = 2 and q — 1. We define 
the partial derivatives of the dependent variables as new variables (prolon- 
gation) in sufficient number to write the equation as second order equation, 
then we can construct a set of differential forms. We speak of the set of 
forms, representing the equations, as an ideal /. It is to be closed. (See [6], 
[7], [11] for more details.) 

In this method Lie derivative of a form with respect to a vector field makes 
an important roll. Lie derivatives of geometrical objects, like differential 
forms, are associated with symmetries of those objects. The Lie derivative of 
a geometrical object carries it along a path, determined by a vector field v, in 
its manifold. If the Lie derivative vanishes, then the vector v represents the 
direction of an infinitesimal symmetry in the manifold. We may construct 
the Lie derivative (£) of differential form in the ideal /. 

Setting the Lie derivative of these forms equal to zero should therefore 
represent symmetries. We require that the Lie derivatives of the forms in / 
to be linear combinations of those forms themselves and when they vanish 
they Lie derivatives also vanish. We can illustrate this by writing C\I = 
mod I, or C V I C /, for Lie derivative properties see [12]. 

1.1 Construction of Differential Forms 

First of all substitute u x by w, thus the KdV equation u xxx + uu x + u t — 0, 
reduced to the following second order partial differential equation 



w xx + uw + u t — 0. 
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The jet space corresponds to the equation (1.1) is a 8— dimensional manifold 
with the coordinate (t, x, u, u t , u x , u tt , u tx , u xx ). Corresponding to the equa- 
tion we have the following contact 1-forms 

9 1 = du — u t dt — u x dx, 

9 2 = du t — u tt dt — u tx dx, 

9 3 = du x — u tx dt — u xx dx. 

Consequently our required forms are 



a 1 
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9 1 A9 2 = (u x u tt — u t u tx )dx Adt + u tx dx A du — u x dx A du t 
+u tt dt A du — u t dt A du t + du A du t , 
a" = 9 1 A 9 3 = (u x u tx — u t u xx )dx Adt + u xx dx A du — u x dx A du x 

+u tx dt A du — u t dt A du x + du A du x , 
a 3 = 9 2 A 9 3 = (u 2 x — u tt u xx )dx Adt + u xx dx A du t — u tx dx A du x 

+u tx dt A du t — u tt dt A du x + du t A du x , 
a A = d9 x = dx A du x + dt A du t , 
a 5 = d9 2 = dx A du tx + dt A du tt , 
a 6 = d9 3 = dx A du xx + dt A du tx , 
a 7 = udt A du — dt A du xx — dx A du. 

Consider a vector field on the assumed jet space in the form of 

, d . d „ d „ d . d d d d 

v = Ai— + A 2 — + A 3 — + A A — + A 5 — + A 6 — + A 7 - — + A 8 - — , 

at ox ou ou t ou x outt ou tx ou xx 

where AiS, i — 1, 8, are smooth functions of (t, x, u, u t , u x , u a , u tx , u xx ). 

The next step is solving the following huge!! partial differential equations 
system, 

CW = Ka\ fori = 1,2, 3, (1.2) 

for smooth functions Aj. 

The system (1.2) has the form 
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d d 
-u(-^—A 2 )u x u tx H (—A 5 ^u tx = 0, 



uu X x UUj tx 

After solving this system with respect to A ± , A 8 , we have 

A 1 = cit + c 2 , A 2 = cix/3 + c 3 t + c 4 , A 3 = -2ci-u/3 + c 3 , 

A 4 = -c 3 u x - dcxUt/3, A 5 = -du x , A 6 = -2c 3 u tx - 8ciu tt /3, 

A 7 = -c 3 u xx - 2c 1 u tx , A s = -AdUxx/3, 

where Ci, c 2 , c 3 , c 4 are arbitrary constants, thus 

/ x <9 ,ci . 9 2ci 9 5ci 9 

V = (dt + C 2 )- + (-X + C 3 t + C 4 )- + (C 3 " T «)^ " + ^U t ) — 

d 8ci d . d Aci d 
-C\U x -= {2c 3 u tx + -z-utt)n [c 3 u xx + 2c 1 u tx )- — u x 



OU x 3 OU tt OU tx 3 OUxx 

We claim that the four vector fields which will be constructed from v, make 
four dimensional symmetry group for KdV equation such as we can obtain 
in Lie method. 

d d 

ox at 

d d d n d d 
v 3 = t— + u x - 2w te - u a 



or aw <7-u t cra ti 

9 9 (9 (9 (9 

v 4 = x— + 3t— - 2u- 5^- 3 u xtt~ 

ox at ou ou t ou x 

9 n d d 
Su tt - Qu t xT: <±u xx - . 

OUtt OUtx OUxx 

First of all we should show that the set of these four vector fields {vi, v 3 , v 3 , v 4 }, 
makes a Lie algebra construction, it is sufficient to show that [vj,v,-] lies in 
the vector space constructed by {v 1? v 2 , v 3 , v 4 }, where [ , ] is the Lie bracket 
of vector fields. By a directly computation we have 

[vi,v 2 ]=0, [vi,v 3 ] = 0, [vi,v 4 ]=vi, 
[v 2 ,v 3 ] = vi, [v 2 ,v 4 ] = 3v 2 , [v 3 ,v 4 ] = -2v 3 , 
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if we label the point parts of these vector fields by 



X - 9 X - 9 

ox at 

d d ^ d „ 9 „ 9 
X 3 = ^tt + 7T' X 4 = x— + 3t- - 2m—, 
ax aw ox at ou 

it is easy to see that the third prolongations of X 1 ,X 2 ,X 3 ,X 4 , vanishes the 
KdV equation, (see [14], [15], for more details.) xf\u xxx + uu x + u t ) = 
0, for i = 1,2,3,4. Where X- , is the third prolongation of X^ thus 
{Xi, X2, X%, X4}, makes a set of four parameter symmetry group for KdV 
equation. It is necessary to say that the four vector fields Vi, V2, V3, V4, which 
we had found are second prolongation of X±, X 2 , X 3 , X±, because in the in- 
cipience of the section, we decreased the order of equation to two. 



2 Solutions and First Integral of KdV Equa- 
tion 

First integral(s) of an ordinary differential equation, which will be defined 
more precisely, is a function which its differential lies in the annihilator of the 
distribution corresponds to the assumed ODE. If we have a first integral(s) 
of an ODE we can find its solution, specially for some differential equations 
which we are unable to find its solution with any model. For the beginning 
of this part, we have some analysis on solutions of equations then, we should 
have a view on some fundamental concepts in theory of differential equations. 
The reader may see [13] for more details. Now consider the KdV equation, let 
G be an n— dimensional Lie group with the Lie algebra g. G is called solvable 
if there exist a sequence of subgroups {e} C G C G 1 C ■ • • C G n _i C G n = 
G, such that each G{ is a normal subgroup of Gj+i. This is equivalence to the 
requirement that the corresponding subalgebras of q satisfies gj+i] C £jj. 



2.1 Bianchi Theorem 

A theorem of Bianchi, [3], states that if an ordinary differential equation 
admits an n— dimensional solvable symmetry group, then its solution can be 
determined, by quadratures, from those to reduces the order of equation; see 
also [15]. 
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Theorem 2.1.1 Let A(x,u^) = be an mth order ordinary differential 
equation. If A admits a solvable n-parameter group of symmetries G such 
that for 1 < i < n the orbits of (ith prolongation of G) have dimension 
i, then the general solution of A can be found by quadratures from the gen- 
eral solution of an (m-n)th order differential equation K{y,w( m ~ n >) = 0. In 
particular, if A admits an m-parameter solvable group of symmetries, then 
the general solution to A can be found by quadratures alone. 

The Bianchi theorem is extended on an arbitrary distribution, which theorem 
(2.1.1) is its special case, [13]. 

By a change of variables as y = x — ct and v = u, the reduced equation 
is v'" + vv' — cv' = 0, where v'" = d 3 v/dy 3 and v' = dv/dy. This can 
be immediately integrated once, v" + v 2 /2 — cv = c\, another integration 
reduced equation to 

111 

-v' 2 + -v 3 - -cv 2 - Cl v - c 2 = 0, (2.3) 
2 6 2' v ; 

where C\ and c 2 are constants. The general solution can be written in terms 
of elliptic function, u = V(x — ct + e), e being an arbitrary phase shift. If 
u — > sufficiently rapidly as \x\ — > oo, then c\ = Ci = in equation (2.3), 
this equation has the solution v — 3c sech 2 \fc/2 + e , provided the wave 
speed c is positive. This produce the celebrated "one soliton" solutions 

u(x, t) = 3csech 2 -y/c(x — ct) + e 

to the KdV equation, which is called Travelling Wave Solution. 

There is a lot kinds of solutions for Kortweg-de Vries, such as Galilean- 
Invariant Solution and Scale-Invariant Solution. 

Consider the flows generated by the four symmetries, 

9i(s)(t, x, u) = (t, x + s, u), 02(s)(t, x, u) — (t + s, x, u), 

9 3 (s)(t,x,u) = (t,x + ts, s + u), 04(s)(t,x,u) = (te 3s , xe s , ue~ 2s ) , 

according to the flows, we can find some special solutions for KdV equation. 
For example if H(x,t) be a solution of KdV, then is so u(t,x) = 5 2 H(5 3 t + 
a, 5x + (3 + jSt) — A, where a, (3, 7, 5 and A are arbitrary constants. We can 
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see that 

u(tx) = ^ t 

u(t,x) = -12tanh(x 2 ) + 8, 

u(t,x) = -127 2 tanh(/5t + 7x + a) 2 + 87 2 --, 
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are another solutions which are constructed with the symmetries and flows. 
If u depends only on x, i,e., u = y(x), then KdV equation reduced to 

y"' + yy' = 0, (2.4) 

and X 1 = d/dx, X 2 = xd/dx — 2yd/ dy, are two symmetries for equation 
(2.4), and x — v(u), y — u, is a change of variable corresponds to X\ for 
equation (2.4). With these new variables and by substituting v' = rj(u) the 
equation (2.4) changes to 

urf + 3r]' 2 - rjrj" = 0, (2.5) 

where rf = dr)/du. This new equation (2.5) has the following three parameter 
symmetry group 

Xi = V 3 -0-, X 2 = r,\^, x 3 = 2u-?--3r ] -?-, 
or] or] ou or] 

corresponding to X 3 , we have the change of variable 

t = ri(u)u*, s(t) = ^\n(u), 

by substituting s' — equation (2.5) reduces to 

4t 4 £ 3 + 3£ - 10t£ 2 + 12t 2 £ 3 + t? = 0, (2.6) 
where £' = d£/dt. This equation has one parameter symmetry group 

X = t(3 + « 2 )^-3(l + t 2 )^- 
The solution of the equation (2.6) is the general solution of KdV equation. 
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2.2 Frobenius Theorem 

The Frobenius theorem is one of the most important theorem in theory of 
differential equations, which one of its results guarantees that an ODE has 
first integral(s). 

Definition 2.2.1 Suppose M is an n— dimensional manifold and p e M. A 
choice of k— dimensional linear subspace D p C T P M is called a k— dimensional 
tangent distribution or a k— dimensional distribution. D p is called smooth if 
D = UpeM D p C TM, is a smooth subbundle ofTM. 

Here T P M and TM are tangent space on M in point p, and TM is the tangent 
bundle on M. 

Lemma 2.2.2 [12] Let M be a smooth n— manifold, and let D C TM be 

a k— dimensional distribution. Then D is smooth if and only if each point 
p G M has a neighborhood U on which there are smooth 1-forms u 1 , ...,u n ~ k 
such that for each q e U , 



then for any uj 1 defined in lemma (2.2.2) we have to 1 E Ann(D). 

If D is a distribution generates by {vi, • • • , v fc }, then D could be discussed 
by {lu 1 , ...,uj n ~ k } too. We will show such a distribution as 



where n — k is codimension of D. 

Definition 2.2.3 Suppose D = J-{vi_, v^) is a k— dimensional distribu- 
tion. The distribution which is generates by the vector fields {vi, Vk] 
and by all possible sorts of commutators [vi, Vj](i < j; i,j = 1, k), is called 
the first derivative of D, i.e., 



D Q = ker u fl • • • D ker u 



,n—k 



More precisely, if we denote the annihilator of D by 




D = ^"(vi, v fc ) 



D (1) = F(v U Vk, [Vi, V 2 ], [V!, V k ], [Vk-l, v k ]). 
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Lemma 2.2.4 [13] Let D = J-{v\, Vk) is a distribution such thatJ r (uj 1 , u n k ). 
Then = D if and only if for i — 1, ...,n — k 

duf 1 A cj 1 A ... A uj n ~ k = 0. 

Definition 2.2.5 A smooth distribution D on a smooth manifold M is called 
Completely Integrable Distribution or a CID, if all points of M contain in an 
integral manifold of D. (A submanifold N C M is called an integral manifold 
ofD ifT p N C D.) 

Now we are ready to give the Frobenius theorem. 

Theorem 2.2.6 [12] Let D is a smooth distribution on a smooth manifold 
M. D is CID if and only if = D. 

Now by using the Frobenius theorem it will be shown that any first order 
ODE has first integral(s). This is the result of the Frobenius theorem and 
its following corollary, before we have a necessary definition. 

Definition 2.2.7 Let D is a smooth distribution on smooth manifold M, a 
smooth function tp G C°°(M) is called a first integral for D if dip G Ann(D). 
Here C°°(M) is the set of all smooth real valued function on manifold M. An- 
other definition for a CID, is that D is CID if and only if there exist n—k first 
functional independent integrals ip ± , p n -k such that D = T{dp\, dp n _k), 

so if a distribution is CID it means that it has first integral(s). 

Corollary 2.2.8 Suppose D is a smooth distribution such that D = JF(cj). 
Then D is CID if and only if 

lu A du = 0. (2.7) 
Now consider a first order ODE 

y' = ^ = f(x,y), (2.8) 

it's clear that the equation (2.8) obtained by taking the 1-form u = dy — 
f(x,y)dx, to zero, so u satisfies the equation (2.7), thus the distribution 
corresponds to the equation (2.8) has first integral(s). 
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2.3 Symmetries of Distribution 



Definition 2.3.1 A symmetry of a distribution is the transformation of the 
manifold M that maps distribution into itself. In other words, a diffeomor- 
phism F : M — > M is a symmetry of a distribution D if F*(D p ) = Dp(p), for 
all p G M. Here F*, is push forward of F. 

Suppose v is an smooth vector field on manifold M, and 9t be its flow, 
then we know 9 t : M — > M induces a diffeomorphism on M. 

Definition 2.3.2 9 t which has defined above is called an infinitesimal sym- 
metry or a symmetry of a distribution D if 9 t along the vector field v consists 
of symmetries of D. i.e, 9*(D P ) = D 9t ^, for all p e M and t. 

Denote Sym(D) the set of all infinitesimal symmetries of distribution D. 

Theorem 2.3.3 [13] Let D be a distribution and X(V) denotes the set of 
all vector fields on D, then the following statements are equivalent. 

i) v e Sym(D). 

u) G X(V) [v, w] G X(V). 

Hi) Vcj G Ann(D) Cyou G Ann(D). 

In this theorem C\tu, denotes the Lie derivative of uj with respect to v. 

Corollary 2.3.4 Sym(D) has a real Lie algebra structure with respect to the 
commutator of the vector fields. 

2.4 Distribution with a Commutative Symmetry Al- 



Let g be a commutative symmetry Lie algebra which its dimension is equal 
to the dimension of J-{uJ l , ...,u n ~ k ) is equal to k. Let {vi, v^} be a basis 
of g and let D is a CID. Then form the matrix 



gebra 




Z 
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because of independence of u/s, then Z _1 is exist, now we are going to con- 
struct a new basis u 1 , ...,u k for Ann(D). These basis are constructed by the 
following relation 



( q1 


\ 




( ^ 


\ 






= z- 1 








) 






) 



It is possible to see that u l s are closed, so the functions tpi(p) = are 
called first integrals of D. Here a is a path from the fix point po to a point 
p G M. Because of closeness of u) l s, the integration is independent from 
choice of a. 

Consider a completely integrable distribution D = T = (ui) with a sym- 
metry v, then according to the relation (2.9), the 1-form Hi = ^fj^, is closed 
and the function 




(2.10) 



is a first integral of D. 

For another example, consider an ODE in the form of equation (2.8), 
we know that the corresponding 1-form is oo = dy — f(x,y)dx, suppose that 
v = a(x, y).d/dx+b(x, y).d/dy, is a symmetry of the above distribution, that 
is the symmetry of the equation (2.8), Z = b(x,y) — f(x,y)a(x,y) and the 
differential 1-form oj = (dy — f(x,y)dx)/(b(x,y) — f(x,y)a(x,y)), is closed, 
and the function 

(p= u, 

J a 

is a first integral of the equation (2.8). The function Z -1 is called an inte- 
grating factor for the equation (2.8). 



2.5 First Integral of KdV Equation 

After we reduced the equation to an ODE, we can obtain its first integral 
due to equation (2.10). The equation (2.3) could be written as 1-form 



ui 



dv — (y\J cv 2 + 2ciV + 2c 2 — v 3 /3jdy, 
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such that its independent variable y, does not enter explicitly, thus it is 
obvious that v = d/dy, is a symmetry for the equation (2.3), consequently 
the closed 1— form cu = —uj / ^ +cv 2 + 2c±v + 2c 2 — v 3 /3, is constructed and 

the function ip = f a (dy — dvj \J+cv 2 + 1c\v + 2c 2 — v 3 /3j, which is a line 
integral on any arbitrary path a is the first integral (dip G Ann(D) where 
D = !F{uj) = J^(v)) for KdV equation, and the set {u(t,x) : p = 0}, is the 
set of all solutions for the equation. 
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